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Abstract 

This is a pilot paper serving as a launching pad for study of orbital and spin 
evolution of binary asteroids. The rate of tidal evolution of asteroidal binaries 
is defined by the dynamical Love numbers ki divided by quality factors Q. 
Common in the literature is the (oftentimes illegitimate) approximation of the 
dynamical Love numbers with their static counterparts. Since the static Love 
numbers are, approximately, proportional to the inverse rigidity, this renders a 
popular fallacy that the tidal evolution rate is determined by the product of 
the rigidity by the quality factor: ki/Q oc In reality, the dynamical 

Love numbers depend on the tidal frequency and all rheological parameters of 
the tidally perturbed body (not just rigidity). We demonstrate that in aster¬ 
oidal binaries the rigidity of their components plays virtually no role in tidal 
friction and tidal lagging, and thereby has almost no influence on the intensity 
of tidal interactions (tidal torques, tidal dissipation, tidally induced changes of 
the orbit). A key quantity that overwhelmingly determines the tidal evolution 
is a product of the effective viscosity t] by the tidal frequency x ■ The func¬ 
tional form of the torque’s dependence on this product depends on who wins 
in the competition between viscosity and self-gravitation. Hence a quantitative 
criterion, to distinguish between two regimes. For higher values of rjx , we get 
ki/Q (X. l/(r?x) ; while for lower values we obtain ki/Q oc rjx- Our study rests 
on an assumption that asteroids can be treated as Maxwell bodies. Applicable 
to rigid rocks at low frequencies, this approximation is used here also for rubble 
piles, due to the lack of a better model. In the future, as we learn more about 
mechanics of granular mixtures in a weak gravity field, we may have to amend 
the tidal theory with other rheological parameters, ones that do not show up in 
the description of viscoelastic bodies. This line of study provides a tool to ex¬ 
ploring the orbital history of asteroidal pairs, as well as of their hnal spin states. 
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1 Dynamics of binary asteroids 

Numerous binary asteroids have been discovered by now (see, e.g., Pravec & Harris 
2007, Taylor & Margot 2011, Jacobson & Scheeres 2011, Margot et ah 2015). Analysis 
of their orbital and rotational dynamics will give us a key to understanding these 
asteroid’s internal properties from observation of their orbits and spin. 

Spin-orbit evolution of a binary is implemented through tidal torques that the com¬ 
panions exert upon one another. Torques due to the permanent figure also participate 
in this interaction and play a role in tidal capture into spin-orbit resonances (e.g., 
Noyelles et al. 2014, Batygin & Morbidelli 2015). In this paper, however, we shall 
concentrate on tides. 

1.1 Difficulties in modeling of binary asteroids’ tidal dynamics 

At hrst glance, the strategy and methods for exploring binary asteroids’ dynamics 
should be the same as those employed in the studies of planets and moons. In reality, 
we are facing a more complicated case, for three main reasons: 

1. Many asteroidal pairs are very close, so calculation of their tidal interaction 
requires inclusion of terms higher than quadrupole (Mathis & Le Poncin-Lahtte 
2009, Compere & Lemaitre 2014, Taylor & Margot 2010). 

2. Tidal reaction of a connected (not rubble) small body is overwhelmingly defined 
by rheology and only to a minuscule degree by self-gravitation (which is not the 
case for planets - see, e.g., Efroimsky 2012b). 

3. Many asteroids, however, are believed to be loose aggregates. Their rheology is 
largely dehned by microgravity, and may differ from the rheology of solids not 
only in parameters’ values. A constitutive equation for granular matter in a weak 
gravity held may, in principle, contain parameters other than rigidity or viscosity. 

1.2 Tidal torques 

When an extended celestial body is deformed by the potential of an exterior perturber, 
the distortion results in an additional tidal potential produced by the extended body 
itself. The latter potential acts back on the perturber and changes its orbital motion. 

In each point of the extended body, the exterior perturbing potential changes in 
time and usually contains an inhnite set of Fourier harmonics. We shall assume the 
reaction linear, so the Fourier spectrum of the reaction (i.e., of the shape changes and 
of the resulting tidal potential) mimics the spectrum of the perturbing potential. 

As ever in physics, the reaction lags behind action. Consequently, one may expect 
that the Fourier modes of the tidal elevation and tidal potential should always lag 
behind the corresponding modes of the perturbing potential. This point, however, 
needs attention; falling behind in time does not £x the sign of a geometric lag. 0 

^ Causality requires reaction to fall behind action, in time. So each Fourier mode of the additional 
tidal potential falls behind the appropriate mode of the perturbing potential, by a time lag At which 
is mode-dependent but is always positive definite. This, however, does not imply that the phase lags 
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Since the perturber “feels” the tides it itself creates in the extended body, the 
perturber experiences an orbital tidal force and an orbital tidal torque. Then the 
deformed body experiences a tidal torque of an opposite sign, which acts to change 
this body’s spin. This interaction is reciprocal, in that the deformed body also acts 
as a perturber upon its companion - and exerts a despinning torque thereon, and 
experiences an opposite orbital torque. 


1.3 A specific problem 


In the literature on tidal dynamics, it has long become a convention to introduce the 
Love numbers and quality factors as separate entities. More often than not, authors 
postulate some frequency-dependence of the quality factors, treating the Love num¬ 
bers as constant parameters (e.g., Mignard 1979, 1980, 1981; Correia & Laskar 2004; 
Efroimsky & Lainey 2007; Cheng et al. 2014). Above that, the dynamical Love num¬ 
bers are very often identified with their static counterparts (Goldreich & Peale 1968, 
Peale & Cassen 1978, Murray & Dermott 1999, Taylor & Margot 2010). 

The static Love numbers are, approximately, proportional to the inverse rigidity. 
For example, the static quadrupole Love number is given by 
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with g, p, and R being the surface gravity, density and radius of a spherical body. 
This altogether renders a common misunderstanding that the tidal evolution rate is 
determined by the product of the rigidity by the quality factor: ki/Q oc. l/{pQ). 

A minor issue with this treatment is that the tidal quality factors, as functions of 
frequency, are different for different values of the degree I (and are different from the 
seismic quality factor). Negligible at ordinary frequencies, these differences become 
critical in the zero-frequency limit, i.e., at crossing of spin-orbit resonances (Makarov 
et al. 2012, Noyelles et al. 2014). 

A bigger issue - the one to be dealt with - is that the dynamical Love numbers, aside 
from their frequency dependence, depend on all rheological parameters of the tidally 
perturbed body (not just rigidity). □ Our goal is to demonstrate that in asteroidal 
binaries the rigidity of their components plays virtually no role in tidal friction and tidal 
lagging. Consequently, it has virtually no influence on tidal despinning, or tidal heating, 
or tidally induced orbital evolution. As we shall see, a parameter that overwhelmingly 
determines the tidal torques (and therefore dissipation, as well as orbit changes) is a 
product of the effective viscosity p by the tidal frequency y . 


should always be positive. Indeed, the phase lag at a certain Fourier mode is a product of this mode by 
the appropriate time lag (see the equation [75] below), so the sign of this lag always coincides with that 
of the mode. Likewise, the observable geometric lags (related to the phase lag through the formula |79| 
below), share the sign with the Fourier mode. Accordingly, no matter whether the perturber is below 
or above the synchronous orbit, some components of the tide, geometrically, fall behind the perturber; 
and some lead. For example, the principal (semidiurnal) Phobos-raised bulge in Mars falls behind. 
Of the lesser bulges, some fall behind and some lead, at various angular rates. 

^ The Andrade model serves as an example of rheology containing parameters other than rigidity 
and viscosity - see the equations (EIH) in Section 13.21 below. More examples of the kind can be found, 
e.g., in Henning et al. (2009). 
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Table 1. Symbol key 


Notation 


Description 
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.the rotation angle of the body 

.the radius of the body 

... .the polar component of the tidal torque acting on the body 
the instantaneous distance between the body and the perturber 

.the semimajor axis 

.the orbital eccentricity 

.the orbital inclination (obliquity) 

.the longitude of the node 

.the argument of the pericentre 

.the mean anomaly 

.the mean motion 

.Newton’s gravitational constant 

.a tidal Fourier mode 

.a forcing tidal frequency 

.a shorter notation for 

.the viscoelastic characteristic time (Maxwell time) 

.the inelastic characteristic time (Andrade time) 

.the viscosity 

.the unrelaxed rigidity modulus 

.the unrelaxed compliance 

.the Andrade parameter 


2 Notation and formalism 

Referring the reader to Appendix for more detail, here we provide basic formulae 
to be used. In Table 1 and thereafter, the word body will mean: the tidally perturbed 
body. Its companion will be termed as the perturber. Evidently, their roles are always 
interchangeable, and the overall dynamics is dehned by tides in both companions. 

2.1 Modes and frequencies 

Both the perturbing potential and the response (tidal distortion of the perturbed body 
and the resulting additional tidal potential) can be expanded over the Fourier modes 

oJijnpq = {I - 2p) u + {I - 2p + q) n + m {(l - 9) ^ {I - 2p + q) n - m 9 (2) 

parameterised with four integers I, m, p, q (Kaula 1964). Here 9 and 9 are the 
rotation angle and rotation rate of the tidally perturbed body, while u and D are 
the perturber’s argument of the pericentre and the longitude of the node, as seen from 
the perturbed body. The “anomalistic” mean motion n is defined through n = M. 
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(with = 0 in the pericentre). The modes can assume either sign, while the 

physical forcing frequencies in the distorted body are positive dehnite: 

Xlmpq = I ^Impq | ~ | (^ “ 2 p + g) Tl - Trt 0 | , ( 3 ) 

SO 

^Impq = Xlmpq ^^^i^lmpq) ' ( 4 ) 


2.2 The despinning tidal torque 


Development of the polar torque can be found, e.g., in Efroimsky (2012 a). The torque 
contains both an oscillating and a secular part, the latter looking as 

(r™) = 
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Here G is Newton’s gravity constant; M is the mass of the perturber; a, i, e are the 
semimajor axis, inclination, and eccentricity; and the angular brackets (...) denote 
time averaging. The notations Fimp{i) and Gipq{e) signify the inclination functions 
and the eccentricity polynomials. The degree-/ Love numbers ki{ujimpq) and phase lags 
(■li^impq) are functions of the Fourier tidal modes ([2]) - see Appendix [Al 

In the series ([5]), some terms are negative (despinning) and some are positive 
(spinning-up). Usually, the leading term is despinning and so is the overall torque. 
This, however, is not obligatory. For example, when the initial spin is prograde but 
slower than the mean motion, or when the initial spin is retrograde, the tidal torque 
acts to increase the spin rate in the prograde sense (Noyelles et al. 2014). Also mind 
that capture into spin-orbit resonances takes place due to a combined action of the tidal 
torque and a torque caused by the permanent triaxiality. Entrapment into higher than 
1;1 spin-orbit resonances is more probable for tighter binaries and higher eccentricities. 


2.3 The orbital tidal torque 

The orbital tidal torque acting on the perturber is opposite to the torque ([5]) wherewith 
the perturber despins or spins up the tidally distorted body. This indicates that the 
products h/Qi = sin €^( 0 ;,^^^) emerging in (0 must be key players not only 

in the spin evolution but also in orbital evolution. 

In practice, orbital calculations employ not the orbital torque, but the tidal poten¬ 
tial. As distinct from the torque, the potential carries the factors ki{uji^^^) cos . 

Recall, however, that the potential is inserted into the Lagrange- or Delaunay-type 
planetary equations. There it is differentiated with respect to orbital elements, and we 
once again end up with[^ ^z^^pg) ■ 

Although our paper addresses hjQi = /^((cu^^pq) sine;(ca^^p^), for completeness we 
also provide in Appendix I A. 6 1 an expression for ki{uji^^^) cos € 1 ( 001 ^^^), for the Maxwell 
body and for some other rheologies. 

In the planetary equations, we also get small terms with hi^impq) cosedw^^p^). They emerge 
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3 Formalism to describe rheology 


As explained in Appendix in an Impq term of the torque we can switch from the 
Fourier mode to the physical frequency Ximpq = I ^impq I ^ 

f^ii^impq) = k^ixi^pg) sine;(x^^pg) . (6) 

The functional form of the dependence of (y) sin (x) upon the frequency y = Ximpq 
is dehned by self-gravitation of the body and by its rheological properties. 

It is common in the literature to denote sin with 1/Qi, and to term the so- 
defined quantities Qi as the tidal quality factors: 


—7--- = sin|e;(a;^^p^)| =smei{ximpq) 

^l\Mmpq) 


(7) 


In this notation, 

^Mmpq) ^^^^Mmpq) 





( 8 ) 


3.1 Constitutive equation 

Rheological properties of a material are encoded in a constitutive equation intercon¬ 
necting the present-time deviatoric strain tensor u^y{t) with the values that have been 
assumed by the deviatoric stress over the time t' < t . Under linear deforma¬ 

tion, the equation has the form of convolution, in the time domain: 

2u^y{t) = J{t) = f j {t -t') a^y{t')dtt' , (9) 

J —00 

and the form of product, in the frequency domain: 

2 = J{x) • ( 10 ) 

Here u^y{x) d'-yy(x) are the Fourier images of strain and stress, while the complex 
compliance J(x) is a Fourier image of the kernel J{t — t') of the integral operator 
(jHl). See, e.g., Efroimsky (2012a,b) for details. 


3.2 Solids. The Maxwell uiodel aud its geueralisatiou 

At low frequencies, deformation of most solids is viscoelastic and obeys the Maxwell 
model. This model can be represented with a viscous damper and an elastic spring 

because the tidally deformed body is no longer a perfect sphere, and its mass cannot be assumed to be 

concentrated at its centre. As a result of this, the anomalistic mean motion n =A4 slightly exceeds 
the Keplerian mean motion G {M + m)/a^ . The terms with ki ) cos e/ ) do not show 

up in the planetary equations for a, e , i, and cause only small secular changes in At, w and . 
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connected in series. Experiencing the same force, these elements have their elonga¬ 
tions snmmed np. This illustrates a situation where a deviatoric stress § generates a 

(e) (v) 

deviatoric strain comprising a purely elastic part U and a purely viscous part U : 


(e) (v) (e) (e) (v) 9 , , 

U = u + U , where § = 2/i U and S = ^r] — U , (11) 


T] and p denoting the viscosity and unrelaxed rigidity. 

As in the Maxwell regime both parts of the strain are generated by the same stress 


{v) (e 

= S = § 


the formulae flTTD can be assembled into one equation: 


or, equivalently: 


U= 7^ § +7^§ 
2/i 2ri 


S + —S = 2/i U , 


with the Maxwell time introduced as 


r,. 


1 


( 12 ) 


(13a) 


(13b) 


(14) 


Comparing the equation flT^ with the general expression ([9]) for the compliance oper¬ 
ator, we see that the kernel is a time derivative of the compliance function 


(Maxwell) 

J{t - t') 


J -\- {t 



t') 


(15) 


0(t — t') being the Heaviside step function, and J being the unrelaxed compliance: 


J = 


1 

h 


(16) 


In the frequency domain, the equation (IT^ can be cast into the form (ITOD . with the 
complex compliance given by 


(Maxwell) — 

J (x) 



J 1 



(17) 


the terms J and — f/(hx) giving the elastic and viscous parts of deformation, corre¬ 
spondingly. Such a body becomes elastic at high and viscous at low frequencies. 

It is often convenient to employ the dimensionless versions of the compliance: 


(Maxwell) — (Maxwell) _ 

J{x)= = l- ^{XTJ-\ 


(18) 
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(19a) 


and of its real and imaginary parts: 

(Maxwell) 

^ix) = 1 


(Maxwell) 

%X) = -(Xrj-^ ■ (19b) 

The Maxwell rheology is a special case of the Andrade rheology that comprises inputs 
from elasticity, viscosity, and inelastic processes (mainly, dislocation unjamming): 

J{X) = J + /3{tx)-^r{l + a) - — (20a) 


= J + /3{tx)-‘^T{l + a) - , (20b) 

where T is the Gamma function, while a and (3 denote the dimensionless and dimen¬ 
sional Andrade parameters. 

While conventional, the expressions fl20al - I20bj) for the Andrade compliance suffer 
an inconvenient feature, the fractional dimensions of the parameter (3 . It was therefore 
suggested in Efroimsky (2012 a, b) to shape the compliance into a more usable form 

Hx) = [l + (*X^A)““r(l + a) - i {x'TmY^] ) (20c) 

with the parameter introduced via 

p = J r;“ . (21) 

and named as the Andrade time. The compliance (I20cl) is identical to (I20ap and (]20bll . 
but is spared of the parameter of fractional dimensions. 

In the dimensionless form, the Andrade complex compliance reads as: 


J (x) = = 1 + {ixr^) “ r(l + a) - i (yr^) ^ 


J 

its real and complex parts being 

S(x) = Tle[J(x)\ = 


( 22 ) 


ne\J(x)] 

J 


= 1 + (xtJ “ cos(^^) r(Q; + 1 ) , (23) 

3(x) - I»P(X)1 = 

= -(XT„)"‘ - sin(^^j r(a + 1) . (24) 

At frequencies below some threshold (Karato & Spetzler 1990, Eqn. 17), dislocation 
unjamming becomes less efficient a process, and the mantle response becomes purely 
viscoelastic. This means that the Andrade time rapidly increases (so the Andrade 
term in the formulae fl20ll falls off rapidly) when the frequency y = Ximpq Soes below 
the threshold value. As a result of this, at low frequencies the complex compliance of 
solids approaches that of a Maxwell body, and we obtain the formulae flTTI- IT^ . 

Mind that in this section the rigidity and compliance assume their unrelaxed values: 

jj, = /i(0) and J = J(0) . (25) 
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3.3 Rubble piles 

Ample observational data on rotation of asteroids have led to a consensus that many 
(possibly, most) asteroids are loose aggregates kept together by gravity forces. Rheol¬ 
ogy of such objects should be borrowed from the recent studies on mechanics of granular 
matter. Such studies render essentially nonlinear constitutive equations which contain 
parameters other than rigidity or viscosity. We however expect that averaging of such 
equations over a volume much larger than a typical size of a granule should furnish 
linearisable equations with effective rigidity and viscosity and, possibly, other effective 
rheological parameters. 0 We leave this part of work for future. 

In this paper, our goal is simply to demonstrate that tidal evolution of an aster- 
oidal binary is parameterised not only by its components’ rigidity, but also by other 
rheological parameters - like the viscosity. 

Partial justification of treating rubble as a viscoelastic material comes from numeri¬ 
cal modeling - see, e.g., the work by Walsh et al. (2008) who studied both monodisperse 
(same-size) piles and simple bimodal distributions (those containing two different sizes 
of particles). It was found that monodisperse aggregates behave in a manner distant 
from fluid, while mixing of spheres of different sizes yields a closer-to-fluid behaviour. 
Within a somewhat different model by Tanga et al. (2009), even monodisperse piles 
demonstrated hydrodynamical behaviour. At the same time, simulated rubble piles 
did not behave like perfect fluids, due to their ability to sustain shear stress. (Also see 
the experimental study by Murdoch et al. 2013.) This means that the aggregates have 
a hnite shear rigidity /i and may be regarded viscoelastic. 


4 Why we need complex Love numbers, 
not just ki and Q 

It has long become common in the astronomical literature to introduce, separately, 
a degree-/ Love number ki and the corresponding quality factor Q. A minor flaw 
in this tradition is that, fundamentally, the inverse of Q is the sine of the degree-/ 
phase lag: 1/Q = sine/, so we better equip it with an appropriate subscript; Qi. A 
bigger inconsistency is that ki and Qi cannot be introduced independently, because 
they both are determined by self-gravitation and rheology of the body. As we shall see 
below, this circumstance becomes crucial in description of spin dynamics of asteroids. 
The only path to consistent treatment lies in introducing complex Love numbers. 

Even if we are aware of the existence and relevance of some rheological parameter, its exact 
physical meaning is inseparable from context, i.e., from a particular way how the parameter enters the 
constitutive equation. For example, the same couple of parameters, rigidity and viscosity, enter the 
Maxwell and Kelvin-Voigt models in different ways. Accordingly, the physical meaning of viscosity 
within these two models is different. Likewise, different is the meaning of rigidity in these models. 
Enough to say, the viscous limit is approached through /r —>■ oo, in the Maxwell model, and through 
/i ^ 0 , in the Kelvin-Voigt model. 

In their often cited elegant work, Goldreich and Sari (2009) calculated the effective rigidity of a granular 
aggregate. Although a big step forward, that project still requires continuation. It is necessary to 
carry out a similar estimate for the effective viscosity of such a medium. An even bigger challenge will 
be to study how exactly these two (and, possibly, other) parameters enter the constitutive equation. 
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4.1 Static tide. Love numbers 


Let a homogeneous spherical incompressible primary of radius R be perturbed by a 
point-like exterior stationary secondary called perturber. Introduce a spherical coordi¬ 
nate system associated with the primary. The latitudes 0, (p* are reckoned from the 
primary’s equator, while the longitudes A, A* are reckoned from a hxed meridian. The 
perturber of mass M is stationarily located in an exterior point r * = {r*, X*, (p*) . 

Consider an arbitrary exterior point r = (r, A, 0) and the surface point R = 
(i?, A, 0) located beneath it (so R < r, r*). In the surface point R, the disturbing 
potential due to the perturber located in r * can be written as 

OO 

W{R,r*) = ^lT^(i^, r*) , with Wi{R, r *) oc Ricosj) , (26) 

1=2 


Pi{cos'y) being the Legendre polynomials, and 7 being the angle between the vectors 
r* and R originating from the primary’s centre. 

Likewise, the additional tidal potential U{r, r*) in the exterior point r, due to 
the primary’s deformation caused by the perturber being in r * , can be expanded as: 

00 

U{r,r*) = ^ Ui{r, r *) , with Ui{r, r *) oc Pi{cos'y) . (27) 

1=2 

A degree-/ spherical component Ui{r, r*) in the point r is proportional to the 
degree-/ spherical component Wi{R, r*) of the perturbing potential in R : 

Ui{r,r*)= kiWi{R,r*) . (28) 


Then the additional exterior tidal potential of the primary can be written down as: 


U{r, r*) 


^ f/i(r, r*) 
1=2 



ki W,{R, r-) 


(29) 


The quantities ki are the static Love numbers. Solution of a boundary-value problem 
for potentials, deformation, and stress renders the following expression: 


k 


(static) 

I 


3 1 

2 (/ - 1 ) ^ ^{static) 


3 1 _ 3 1 

2(/ - 1) 1 + Bifi ~ 2(/ - 1) 1 + Bi/J 


where 

(2/2 + 4 / + 3) 3 (2/2 + 4 / + 3) 

IgpR ~ AlnGp^R^ 


(31) 


and 


(static) 

Ai = Bip 


Bi/J 


(32) 


^ Since the points r and R are located on the same ray pointing out of the centre of the body, 
the angle between r * and r is equal to the angle 7 between r * and R. 


(30) 
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In the above formulae, p, g, and R are the density, surface gravity, and radius of the 
body, while G signifies the gravitational constant. Since static implies: relaxed, then 
the rigidity modulus p and the compliance J = 1/p assume their relaxed values: 


p = p(oo) and J = J(oo) 


in the static case only. 


(33) 


This will not be the case for evolving tide to be considered below. There, p and J 
will assume their unrelaxed values - like in Section 13.21 

To draw this section to a close, we would mention that in the denominator of (l30l) 

(static) 

the term 1 comes from self-gravitation, while the term Ai originates owing to 

(static) 

rheology. This way, Ai serves as a dimensionless measure of strength-dominated 
versus gravity-dominated static behaviour; see Section 15.2.11 below. 


4.2 Evolving tide. Love operators and complex Love numbers 

In the case of a perturber “feeling” the tide it itself creates, we set r = r*. This 
enables us to use a shorter notation U{r) instead of U{r, r*). When the tide is 
evolving, the time-dependence has to be introduced, both in the overall tidal potential 
U{r, t) and in its Legendre components Ui{r, t). 

Under evolving tidal disturbance, the expression fl29|) becomes a linear operator 
(called the Love operator): 


' R \ R 


(34) 


where the Love functions ki{t — t') are analogues to the compliance functions J{t — t') 
introduced in the equation (j9]). Overdot stands for time derivative. 

In the frequency domain, the convolution fl34l) becomes [§ 


_ / R \ _ 

C(,(X') = (7] k(x)Vi(x) 


(35) 


with X being the frequency. Here (y) and H/ (y) are the Fourier or Laplace com¬ 
ponents of the potentials U; it) and Wi (t ), while the complex Love numbers (y) are 
the Fourier or Laplace components of the kernels (t — f'). 

The frequency-dependencies fc;(y) should be derived from the expression for the 
complex compliance </(y) or from its inverse, the complex rigidity p(y) = 1 /J(y). 
Into that derivation, also the size and mass of the body will enter. So the tidal response 
will be defined by both rheology and gravity. 

Referring the reader to Efroimsky (2012a,b) and the literature cited therein, we 
state that for viscoelastic bodies the complex Love number ki{x) is related to the 


® The equation (1351) is simplified, in that a more accurate treatment would render the spectral 
components as functions of the tidal mode uj = and not of the physical frequency y = Ximva ~ 

l^impq I • Fortunately, though, the simplification is affordable. As explained in Appendix IA.4I below, 
the employment of the positive definite physical frequencies instead of the tidal modes is legitimate, 
insofar as we do not forget to include the sign multiplier Sgnwj^p^ into an Impq term of the expansion 
for the torque. With this stipulation in mind, in this article we shall always expand over y. 
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complex rigidity fi{x) by the same algebraic expression through which the static Love 
number ki is related to the relaxed rigidity p(oo). Equivalently, one can say that the 
complex Love number ki{x) is related to the complex compliance J(x) = l/fl{x) by 
the same algebraic expression through which the static Love number ki is related to 
the relaxed compliance J{oo) = l//i(cxo) : 0 


kiix) = 


1 + Bikiix) 2{l-l) 1 + Bi/J{x) 

Writing the degree-/ complex Love number as 

h{x) = + ilm[ki(x)] = |fc,()()| e 


(36) 


(37) 


we conventionally denote the phase as — , with a “minus” sign. This convention 

imparts ei with the meaning of phase lag, as can be understood from the relation ([35]). 
Equating the absolute values of the right- and left-hand sides of the relation (l35i) . we 
make another useful observation; the role of the dynamical Love number is played by 
the absolute value of the complex Love number: 


h = \h(x)\ 


(38) 


Finally, we notice that 


\ki{x)\ sinei(x) = -Xm[ki{x)] 

whence, using the expression (j36|) . we arrive at 


(39) 


ki{x) sinez(x) 


3 _ B, Xm [J{x)] _ 

2(/ - 1) (7^e [J{x)] + B,Y + (Xm [J{x)]f 


(40a) 


where ki{x) = l^((x)l cind, as ever, y is a shortened notation for Ximpq ■ 

We see that the product ki{Ximpg)/QiiXimpq) = kiiXi^pg) smciiximpg) standing 
in an Impq term of the expansion for the tidal torque can be calculated consistently 
through the formulae fl5B] - l40a]) . with the coefficients Bi rendered by fIM]) . For a homo¬ 
geneous incompressible sphere, the only information needed to calculate this product 
is: the radius, the density, and the rheological law J(y). 

The above expression can be cast into the form of 


ki{x) sine;(y) = 


where 


3 _ A ^(x) _ 

2(^-1) (3?(y)+A)' + »(x)^ 


(40b) 


^{x)^ne[J{x)]/J , ^(y)^Xm[J(y)]/J , (41) 

^ This algebraic similarity is but one manifestation of a more general theorem known as the 
correspondence principle or the elastic-vis coelastic analogy. Although the discovery of this principle 
is sometimes attributed to Biot (1954), an early version thereof can be found in the cornerstone work 
by Sir George Darwin (1879). 
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are the real and imaginary components of the dimensionless complex compliance, while 


A, = BijJ 


(42) 


is the dimensionless rigidity. The dimensionless rigidity Ai relates to Bi in the same 
way as it used to do in the static case, equation fl52D . Mind, though, that now p and 
J denote the unrelaxed values of the rigidity and compliance, not the relaxed ones as 

(static) 

in the static problem. So the values of Ai are different from Ai 

While the product /QiiXimpq) = ^liXimpq) ^^^^liXimpq) shows up in an 

Impq term of the tidal torque, a product 


kiiXin^pq) coseiiximpq) 


IZe 



(43) 


emerges in an Impq term of the expansion for the tidal potential. For reference pur¬ 
poses, we provide in Appendix IA.6l a general expression for this product, as well as its 
versions for the Maxwell and other rheologies. 

Do draw this section to a close, it would not hurt to mention that the expression 
fl36|) . as well as the ensuing formulae fHOj) . are inapplicable to fluid celestial bodies.!^ 


4.3 The tidal lag versus the “seismic” lag 

From the expression flH^ . we obtain the actual dynamical Love numbers: 


ki{x) = I ki{x) I = 




(»(x)+A,f+Hxr 

The ratio of the formulae fl40bp and fl44p gives us the inverse tidal quality factors: 

A Ax) 


(44) 


Qiix) = sinez(x) = - 


AAx)^ + Ax)^ \/(Ax) + A)^ + Ax)^ 


(45) 


The tidal quality factors should not be confused with the “seismic” quality factor 
whose inverse is the sine of the phase lag 6 showing up in the rheological law: 


Qix) = sin(5(x) = - 


^(x) 


The former and latter formulae interrelate the tidal and “seismic” phase lags: 

.4, 


(46) 


sine/(x) = 


\/( Ax) + + A(x)" 


sin(5(x) . 


(47) 


To the best of our knowledge, this relation has never appeared in the literature hitherto. 

® Tidal response of fluid objects (like Jupiters or stars) cannot be treated as a viscous limit of 
the Maxwell rheology. Dissipation in these bodies is ruled by the equations of hydrodynamics, with 
turbulence taken into account (Ogilvie & Lin 2004; Remus et al. 2012; Auclair-Desrotour et al. 2015). 
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As one would expect, the lags coincide when rheology “beats” gravity - i.e., when 
Ai 2 > I 3?(x) I , I Q'(x) I • Otherwise, the tidal phase lags deviate considerably from 
the nominal “seismic” lag 5. This happens on approach to a spin-orbit resonance, 
i.e., to a zero frequency. While the “seismic” lag stays constant (and approaches ti/2 , 
in the case of Maxwell body), the tidal lags continuously go through zero - and so 
do the products sine^, as can be seen from Figure [H In this hgure, we depicted 
ki sin not as a function of the positive dehnite physical frequency y = Ximpq ’ 
as a function of the Fourier mode cu = that can have either sign.l^ The shape 

of this function is about the same for all linear rheologies. For a Maxwell body, the 
extrema are located at 


^peak 


± 


1 


^ SirGp^ 
57 T] 


(48) 


The peaks’ amplitude is virtually insensitive to our choice of the viscosity rj. At the 



- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 

tidal mode co 


Figure 1 A typical shape of the quality function ki{uj) sine)(a;), where to is a shortened 
notation for the tidal Fourier mode . (From Noyelles et al. 2014.) 

same time, the spread between the extrema depends on 77 . As we see from the formula 
(H 8 l) . for a higher viscosity the peaks are located close to the point of resonance. As the 
viscosity assumes smaller values, the peak frequency spreads out, eventually reaching 
and bypassing the orbital frequency n . In realistic settings, however, this would require 
extremely low viscosities. 

Outside the inter-peak interval, the function ki{u}) sine)(a;) behaves as the inverse 
Q in seismology, and changes very slowly with frequency. 

Between the peaks, the products ki{u)) sine;(a;) are almost linear in cu = , 

and go continuously through zero. This enables the tidal torque to transcend spin-orbit 
resonances continuously (Makarov & Efroimsky 2013, Makarov et al. 2012, Noyelles et 
al. 2014). 

It can be demonstrated that the linearity of ki sine/ in u is equivalent to the 
time lag being frequency independent: Ati{ui^p^) = At, see Efroimsky & Makarov 

9 Recall that Xi„,pg = I ^impg I ■ 
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(2013). 0 This means that the tidal response of an actual terrestrial body can be 
approximated with the constant-At model only when all considered tidal frequencies 
are lower than \ Upeak \ ^ or, roughly speaking, when all mean motions and spin rates 
are lower than | Upeak \ ■ To achieve this, we would have to endow the viscosity with very 
low values. This observation absolutely prohibits the application of the CTL (constant 
time lag) tidal model to solid or semi-molten silicate planets (while for liquihed planets 
this entire formalism is not intended anyway). 


5 Complex Love numbers for different rheologies 

For homogeneous near-spherical objects, the formulae fHOj) provide a turn-of-the-crank 
method. Inserting therein the rheology J{x), along with the values of the radius and 
density, we obtain the functions ki{x)/Qi{x) = k^i^x) sine,(x) , with y = Ximpq- 
The complex compliance J{x) is a function of the tidal frequency and rheological 
parameters. This is why the tidal torque is sensitive both to the frequency and the 
rheological parameters (and, thereby, also to the temperature). A similar conclusion is 
valid for the tidal heating rate (Henning et ah 2009, Makarov & Efroimsky 2014). 


5.1 Self-gravitating Maxwell sphere 


Insertion of the complex compliance flTTll into the formulae fl36l) and 

\{X) = ^ ^ 




1 + * X 

1 


2 (/-l) 1 ^ ^ 


1 + * X 


and 


^ =inqx) = 


1 

‘ vx 


2(1-1) (J + B,)^ + ^ 

'I X 




1 


Fm X 


2(/-l) (1 + AiY + 


M A 


furnishes: 

(49a) 

(49b) 


(50a) 


(50b) 


Given by fl^ . the factors Bi are dimensional, while Ai = Bi ^ = BJ J are dimen¬ 
sionless. Recall that, in distinction from the static case, here jj, and J = l/p denote 
the unrelaxed (not static) values of the rigidity and compliance, wherefore our Ai 

{static) _ 

should not be misidentified with their counterparts Ai from Section 14.1[ 

Generalisation of the formula fl50|) to the Andrade model can be found in Efroimsky 
(2012a, eqn. 170). 

Some authors call the constant-At tidal model viscous. This is not a good term, since in the 
viscous limit of the Maxwell model we still arrive at a kink shape of kfuj) sined^)) like in Figure [TJ 
The CTL (constant time lag) tidal model is a better term for this approach. 
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5.2 Three extreme cases 

5.2.1 Self-gravitating drop of an ideal fluid. No rheology at all 

In the “no-rheology” limit, the material becomes an ideal fluid wherein shear deforma¬ 


tion entails no stress. This implies: 

77 —)-0 , /i—)-0 , J = - >00 , (51a) 

h 

wherefrom J{x) cxo and .4,^ —?■ 0 , (51b) 

though the real part of the dimensionless compliance stays flnite: 3?(x) 1. Thence, 

^ ^( 7 ^ . ( 52 ) 

klix) sine;(x) 0 . (53) 


From fl52D we see that in this limit the complex Love number becomes real. Due to 
the absence of rheology as such, this value of the Love number originates solely from 
self-gravitation. From fl53l) we understand that all terms in the expression for torque 
become zero. No heat is produced either. 


5.2.2 Self-gravitating elastic body 

This limit is implemented through 

^ 00 , r] ^ 00 , with J staying flnite . (54a) 

Thence the complex compliance becomes real: 

J(X) ^ J . (54b) 


This again renders the complex Love numbers real: 


kiix) 


3 1 

2 (/ - 1 ) I + Bifx 


3 1 

2(/-i) 1 + A 


^lix) sine,(x) -)■ 0 . 

Phase lags vanish, no tidal torque emerges, and no dissipation is taking place. 


(55) 

(56) 


5.2.3 Self-gravitating perfectly viscous body 

In this situation, we have to get rid of rigidity by setting /i —)■ cx) (yes, 00 ): 
) p—j-cx) , J = - ?-0 , 

/X 

— % 

wherefrom 4 (y) —)■ — — and A, ^ 00 . 

’K ‘ 


(57) 

(58) 
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( 59 ) 


In this limit, the expressions fl49ap and fISUap become: 


kiix) 


2 (/ - 1 ) 1 + 17^ xBi 


kiix) sine,(x) ^ 


B — 

I VX 


BiVX 


2(^ - 1) R2 + ^ 


2 (/-l) + 1 


(60) 


It should however be noted that the viscous limit in these formulae cannot be inter¬ 
preted too literally. Recall that all our algebra rests on the cornerstone fact that the 
complex ki{x) relates to the complex J{x) in the same manner as their static coun¬ 
terparts do. This is an approximation that stays valid insofar as the accelerations and 
inertial forces can be neglected in the equation of motion (see Efroimsky 2012 a and 
references therein). In the absence of rigidity, the Coriolis force comes into the picture 
and entails complicated hydrodynamics (Ogilvie & Lin 2004, Auclair-Desrotour et ah 
2015). So the compliance J cannot be set exactly zero. It should sooner be understood 
that the formulae fl59l) and fl60l) pertain to a finite J r]x ■ If such highly viscous 
objects exist in the world of small bodies, then the tidal torque (and dissipation) for 
these objects is proportional not to the rigidity fi but to 7]x , as can be seen from the 
above expression. 


6 Numbers 

Given our sparse knowledge of the actual parameters and structure of asteroids and 
comets, we are not sure how good approximation the Maxwell model is. We shall 
employ it for the time being, due to the lack of anything better. 

As for the parameters’ values, we only know that the mean density of most asteroids 
is about p = 2 X 10^ kg/m^ , up to a factor of 1.5 (Pravec & Harris 2007), while the 
density of comets is around 0.6 x 10^ kg/m^ (Britt et ah 2006). 

6.1 Rigidity beats gravity — and goes out 

Assuming the nominal values ^^(cxo) = 5 X 10^° Pa and p = 3.5 x 10^ kg/m^, 
Goldreich & Sari (2009, Section 3.1) obtain, for rubble piles: 


(static) 

A 2 oc 

10 ^ 

for 

R 

10 "^ m , 

(Ola) 

(static) 

A 2 oc 

10 ^ 

for 

R ^ 

10 ® m . 

(61b) 


Although the unrelaxed rigidity differs from its relaxed (static) counterpart, for regular 
solids the difference is not dramatic. We presume that for rubble piles, too, the said 
difference does not exceed two orders of magnitude. If this presumption is right, we 
can accept, for the unrelaxed dimensionless rigidity, that 

> 1 , for R < 10® m . (62) 
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This should not come as a surprise, because Ai being larger than unity implies that 
in small bodies “rigidity beats gravity”. 

Now we can safely assume in the formulae fH9] - l50|) that, for all asteroids and 
comets (and for not too large satellites), 1 + Ai ^ Ai and, equivalently, J + Bi ~ Bi. 
Hence the following approximation for the squared Love number: 



3 

2 1 + X^ 

M X ^ 

3 

-b v"^ x^ 

to 

1 

1 tMi + A)2 

to 

1 

1-^ 

1 -b 772 Bj 


and for ki/Qi : 


kiix) sine;(x) 


B, 


VX 


B, 


_x_ 

VX 


2 (/- 1 ) (J + Bif + 


^2 ^2 


2{l - 1) Bf + ^ 


^2 ^2 


(63) 


(64) 


We see that the Love number becomes rigidity-independent for rjx ^ ■ Moreover, 

our second entity, ki/Qi , does not depend upon rigidity at all.^i] 

The latter fact warrants a seemingly paradoxical conclusion that the rigidity of 
small bodies is too large to be relevant for description of tidal lagging in these bodies; 
too large - even with a reduction derived by Goldreich & Sari (2009). In reality, this 
conclusion works also for monoliths, because in the expression for Ai a higher rigidity 
of a solid rock is partially compensated by its higher density squared. So the conclusion 
works for all small bodies, including satellites, - provided they are (a) not too large, 
so the inequalities fl62D hold; and (b) may be regarded as Maxwell bodies. 


6.2 Viscosity versus self-gravitation. A practical criterion 


Our next goal is to examine the denominators of the right-hand sides of the expressions 
fl63ll and fIMl) . In both denominators, it is important whether Bixv is much larger or 
much smaller than unity. With the expression flTIl) for Bi kept in mind, this furnishes 
us with the following criterion. 

For tidal frequencies and effective viscosities obeying 


57 X V 

Stt Gp2i?2 


> 1 


the expressions for ki and ki/Qi read as: 


0 


(65) 


kiix) 


2{l-l) 


yi + V^X^J^ 

Si VX 


( 66 ) 


11 


We can use the equation (f64|l to express ki/Qi as a function of the dimensionless rigidity Ai : 


kiix) sinedx) = 


A, 


TmX 


- » V + ^ 

' M ^ 


However, since Ai = Bi fi and = g/g, casting our expression into the above form will only 
camouflage the absence of an actual dependence upon g ot J . 


Evidently, for even higher viscosities, the Love numbers become simply ki 


1 


2(1-1) ■ 
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klix) sine,(x) 




3 1 


(67) 


They are good everywhere except in the low-frequency limit. There, the inequality 
(16^ is no longer satished, and we have to switch from fl671) to the full expression fl64)) . 
If the product y r] assumes low values, we arrive at a situation with 

57 X V 
Stt Gp^R^ 

so the expressions of our concern become: 



kiix) ~ 


3 

2 (/-l) 


(69) 


^/(x) sine^(x) ^ 


3 

2 (/-l) 


Bivx 


(70) 


The latter expression can be used to explore tidal entrapment into spin-orbit states. 


7 Conclusions 

In the article, we have shown that in asteroidal binaries the intensity of tidal in¬ 
teractions (i.e., the tidal torques and the ensuing spin and orbit evolution) depends 
primarily on the product of the tidal frequency y by the viscosity p. Contrary to a 
long-entrenched opinion, the rigidity p exerts virtually no influence on tidal evolution. 

We have developed a quantitative criterion that enabled us to separate two types of 
the functional dependence of kJQi upon the said product xh • For higher values of this 
parameter, we obtained: ki/Q (x l/{px) > see the equations ([65]-[67]). For lower values 
of XV 1 got: ki/Q oc pXi see the formulae ([68]-[70D. The dependencies fl67|) and 
(170|) determine the spin histories, including tidal captures into spin-orbit resonances. PI 
Our study rests on an assumption that small bodies can be described with the 
Maxwell model. The model is well applicable to rigid rocks at low frequencies, so 
our conclusions are safely valid for monoliths. We apply this approach also to rubble, 
for the lack of a better model so far. As we learn more about mechanics of granular 
mixtures in a weak gravity field, we should be prepared to amend the tidal theory with 
other rheological parameters, ones that do not show up in viscoelastic mechanics. 

It should be finally added that the old misbelief about the allegedly key role of the 
rigidity has pervaded also the literature on inelastic relaxation of processing (wobbling) 
asteroids and comets (Prendergast 1958; Burns & Safronov 1973; Efroimsky 2001, 2002; 
Molina et al. 2003; Ryabova 2004; Sharma et al. 2005; Breiter et ah 2012; Breiter 
& Murawiecka 2015). In an upcoming work (Frouard & Efroimsky 2015), we shall 
demonstrate that in the said setting, too, the leading role belongs to viscosity. 

The numerical factor on the left-hand side of (l68]) corresponds to 1 = 2. If higher degrees are 
included, the factor should be chosen as 3 (2F -|- 4/ -f 3)/ (AItt) for the maximal I taken into account. 

Be mindful that the approximation (1671) cannot be used in the vicinity of a resonance where x 
vanishes. In that vicinity, one has to employ the full expression (1641) or its approximation (I70p . 
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A Appendix. 

In this Appendix, we briefly summarise the basic facts from the theory of bodily tides. 
For more detailed explanation, see Efroimsky & Makarov (2013) and references therein. 

A.l Fourier modes and physical frequencies 

As we mentioned in Section 12.11 both the potential from the perturber and the addi¬ 
tional tidal potential of the perturbed body can be expanded over the Fourier modes 

^impq - {I-2p + q)n - me , (71) 

with n = M. being the “anomalistic” mean motion and 9 being the spin rate of 
the tidally perturbed body. These modes can assume either sign, while the physical 
frequencies of stresses and strains in the perturbed body are positive dehnite: 

Ximpq = I ^Impq I ~ \ {I - 2p + q) u - m 9 \ . (72) 

A.2 Phase lags and dynamical Love numbers 

The tidal reaction always lags in time, compared to action. Each Fourier mode 
has a time lag of its own, . Owing to causality, time lags are positive dehnite: 

> 0 ■ (73) 

For physical reasons, the time lags are dehned by the forcing frequencies Ximpq = 

I ^impq I ! ignorant of the sign of the modes : 

^tli-^lmpq) = = AtliXimpq) ■ ( 74 ) 

In the Darwin-Kaula theory of tides, the phase lags are products of the modes ooimpq 
by the corresponding time lags: 

^li^lmpq) = ^Impq ^^li^lmpq) (75a) 

= Ximpq ^tliXlmpq) ^S^^i^lmpq) ' (75b) 

The functional forms of the lags and Love numbers (as functions of the Fourier mode 
(^impq) are dehned only by the degree I. Hence the notation: /c;(cu,^pg), 6 ,( 0 ;,^^^), 
Ati{uimpq) ■ This, however, is true only for a homogeneous spherical body. 0 

For oblate homogeneous bodies, the functional dependencies of the lags and Love numbers on 
modes are defined by both the degree I and order m (Dehant 1987, Ogilvie 2013). So one should 
write: fcj^K^pg), eimi^impq) ^ and Aq^(w;„pg). These Love numbers and lags differ from the 
Love numbers and lags of the spherical reference body by terms of the order of the flattening, so a 
small non-sphericity can be neglected. 
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Keeping in mind that the time lag is positive definite, we see from fl75ap 

that the signs of e/ and always coincide. So the phase lags are odd functions of 

ujimpq and thus can be written as 

= \^Mmpq)\ i^lmpq) 

or, equivalently, as 

^Mmpq) = ^liXlmpq) ^S^^i^lmpq) > ( 76 b) 

with eiixi^pg) = \ei{ui^pg)\ being positive-definite. 

Just as the time lags, the dynamical Love numbers depend on the actual physical 
frequencies Ximpq ■ This is the same as to say that they are even functions of oJi^^pq ■ 

f^li-^lmpq) = kMmpq) = ^liXl^pq) • (77) 

From the latter two equations, it ensues that 

f^Mmpq) ^^^^Mmpq) = ^liXlmpq) ^liXlmpq) i^lmpq) ' ( 78 ) 


A.3 Observable geometric lags 

Geometrically, tides comprise a system of superimposed bulges moving around the 
body in different directions, at different angular, and with different lags. The principal, 
semidiurnal bulge has an elliptic shape. As seen from the surface, it moves at the same 
angular rate as the perturber, with a positive or negative lag. Lesser bulges have more 
involved shapes and different velocities. By geodetic measurements, it is possible to 
single out those bulges, and to measure their magnitudes and geometric lags (Benjamin 
et ah 2006, Fanes & Bettadpur 1996). 

The geometric lag angle of an Impq bulge is (Efroimsky & Makarov 2013, Eqn 26): 


5 


Impq 




(79) 


For example, the geometric lag of the principal bulge is ^2200 = ^^2 = {n — 

9) At 2 , where the time lag is taken at the appropriate, semidiurnal mode: At 2 = 
At 2 ( 1 ^ 2200 ) • Also mind that the tide components with m = 0 are longitude-independent 
and the notion of geometric lag becomes meaningless - see the equation (25) in Ibid. 


A.4 Why positive frequencies are insufficient 

While the general formula for a Fourier expansion of a field includes integration or 
summation over both positive and negative modes, in the case of real fields it is sufficient 
to expand over positive frequencies only. 0 Surprisingly, the tidal theory makes a rare 

The condition of a field being real requires that the real part of a Fourier term at a negative 
mode is equal to the real part of the term at an opposite, positive mode. Thus one can get rid of the 
terms with negative modes, at the cost of doubling the appropriate terms with positive modes. (The 
convention is that the actual field is the real part of a complex expression.) 
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exception from this rule. This theory cannot be expressed via the physical frequencies 
Ximpq ~ I ^impq I solely, with uo mentioning of the modes ■ It turns out that a 
contribution from a Fourier mode into the additional tidal potential is not physically 
equivalent to the contribution from the mode of an opposite sign. 

To understand the reason for this, recall that the expansion for the tidal potential 
is needed to obtain expressions for tidal force and torque. These are inhnite series 
too; and in each of these two series an Impq term contains explicitly the sign of the 
appropriate tidal mode . Specihcally, an Impq term in the expressions for the 

tidal force or torque contains a factor 

^li^lmpq) = KUlmpq) ^liXl^pq) i^lmpq) > (§ 0 ) 

where we made use of the formulae fl76bp and fl77|) . 

This way, if we choose to expand the tidal torque or the tidal force over the posi¬ 
tive dehnite frequencies Ximpq shall have to insert “by hand” the multipliers 

Sgnw^^pg into each Impq term of these expansions (Efroimsky 2012a,b). 

A.5 A typical shape of the function k^{oJi^pq) sin 

The function sine;(a;;^p^) = k^iximpq)/QiiXimpq) ^?>^i^impq) enters the 

Impq term in the expression ([S]) for the tidal torque. The kink shape of this function is 
amazingly robust: it is about the same for all linear rheological models, including the 
most extreme ones (say, purely viscous mantle with no rigidity). It will exhibit a sharp 
kink, as in Figured] with two peaks having opposite signs. The fact that the function 
goes continuously through zero enables the torque to transcend spin-orbit resonances 
continuously. 

Close to = 0, the frequency dependence makes two extrema, and demon¬ 

strates a slow frequency-dependence at higher frequencies. 


A.6 Expression for the function cos 

As was mentioned in Section above, a product 

kliXlmpq) ^^^^liXlmpq) = (81) 

stands in an Impq term of the expansion for the tidal potential U. It is equal to 

3 {ne [J{x)] +B,)ne [J{x)] + (Xm [J(x)] Y 


kiiXimpq) ^oseiixi^p^) = 


2(/ - 1) {ne [J(x)] + + (Xm [%)] )' 

(5J(x) + AJ 3?(x) + ^5 (x)2 

..o 


2{‘-l) (S(x) + At + 

For a self-gravitating Maxwell sphere, this expression becomes: 


kliXlmpq) ^^^^l^Xlmpq) ~ 


( 


2(/-l) 


1 - 


A; 




V 


Af H—2^ 


(83) 


(82a) 


(82b) 
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Of a special interest are the three extreme cases considered in Section 15.21 


self-gravitating ideal-fluid drop : 

^liXlmpq) 

^O^^liXlmpq) 

-)■ 

3 

2{l- 

1 ) ’ 

(84) 

self-gravitating elastic sphere : 


^O^^liXlmpq) 


3 

2{l- 

A ) > 

(85) 

self-gravitating viscous sphere : 

^liXlmpq) 

^O^^liXlmpq) 

-)■ 

3 

2{l- 

1 

1 ) 1 + 

(86) 


Recall, however, that the Kaula theory of bodily tides becomes inapplicable in the limit 
of a perfectly viscous body, as was explained in Section 15.2.31 
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